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PROBLEM: OPERATOR ERROR ESTIMATES IN L 2 



T. A. SUSLINA 

Abstract. Let O C R d be a bounded domain of class C 2 . In the 
Hilbert space 1/2(0; C"), we consider a matrix elliptic second order dif- 
ferential operator Ad.e with the Dirichlet boundary condition. Here 
e > is the small parameter. The coefficients of the operator are 
periodic and depend on x/e. A sharp order operator error estimate 
H-A^ 1 — (A%,)~ 1 \\l-,-*L; < Ce is obtained. Here A° D is the effective op- 
erator with constant coefficients and with the Dirichlet boundary con- 
dition. 



Introduction 

The paper concerns homogenization theory of periodic differential oper- 
ators (DO's). A broad literature is devoted to homogenization problems in 
the small period limit. First of all, we mention the books [BeLPa], [BaPan], 
[ZhKO]. 

0.1. Operator-theoretic approach to homogenization problems. 

In a series of papers [BSul-3] by M. Sh. Birman and T. A. Suslina a new 
operator-theoretic (spectral) approach to homogenization problems was sug- 
gested and developed. By this approach, the so-called operator error esti- 
mates in homogenization problems for elliptic DO's were obtained. Matrix 
elliptic DO's acting in £2(1^; C n ) and admitting a factorization of the form 
A £ = 6(D)*g(x/e)6(D), e > 0, were studied. Here <?(x) is a periodic matrix- 
valued function and 6(D) is a first order DO. The precise assumptions on 
g(x) and 6(D) are described below in Section 1. 

In [BSul-3], the equation A E u £ + u e = F, where F € L 2 (M d ;C n ), was 
considered. The behavior of the solution u e for small e was studied. The 
solution u e converges in L2(M rf ; C n ) to the solution uo of the "homogenized" 
equation „4°u + u = F, as e — s> 0. Here A = 6(D)* g°b(D) is the effective 
operator with the constant effective matrix g°. In [BSul], it was proved that 



l u e - u o|lL2(R d ) - Ce \\ F 



\L 2 



2000 Mathematics Subject Classification. Primary 35B27. 

Key words and phrases. Periodic differential operators, homogenization, effective op- 
erator, operator error estimates. 

Supported by RFBR (grant no. 11-01-00458-a) and the Program of support of the 
leading scientific schools. 

1 



2 



T. A. SUSLINA 



In operator terms it means that the resolvent (A £ + I) -1 converges in the 
operator norm in L2(R d ;C n ) to the resolvent of the effective operator, as 
e — s- 0, and 

\\(A e + I)- 1 -(A° + I)- 1 \\ L2m ^ L2m <Ce. (0.1) 

In [BSu2], more accurate approximation of the resolvent (A £ + I)^ 1 in 
the operator norm in L2(M. d ; C n ) with an error term 0(e 2 ) was obtained. 

In [BSu3], approximation of the resolvent (A £ + I)^ 1 in the norm of op- 
erators acting from L 2 (R d ; C n ) to the Sobolev space iJ 1 (R d ; C n ) was found: 

\\(Ae + - (A^ + I)- 1 - sK(s)\\ L2{m ^ mm < Ce; (0.2) 

this corresponds to approximation of u £ in the "energy" norm. Here K{e) 
is a corrector. It contains rapidly oscillating factors and so depends on e. 

Estimates (0.1), (0.2) are called the operator error estimates. They are 
order-sharp; the constants in estimates are controlled explicitly in terms of 
the problem data. The method of [BSul-3] is based on the scaling transfor- 
mation, the Floquet-Bloch theory and the analytic perturbation theory. 

0.2. A different approach to operator error estimates in homogeniza- 
tion problems was suggested by V. V. Zhikov. In [Zhl, Zh2, ZhPas, Pas], 
the scalar elliptic operator — div g(x/e)V (where g(x) is a matrix with real 
entries) and the system of elasticity theory were studied. Estimates of the 
form (0.1), (0.2) for the corresponding problems in M. d were obtained. The 
method was based on analysis of the first order approximation to the solution 
and introducing of an additional parameter. Besides the problems in M. d , 
homogenization problems in a bounded domain O C M d with the Dirichlet 
or Neumann boundary condition were studied. Approximation of the solu- 
tion in H 1 ^) was deduced from the corresponding result in M. d . Due to the 
"boundary layer" influence, estimates in a bounded domain become worse 
and the error term is 0(s l l 2 ). The estimate ||u e — Uo||l 2 (o) — Ce 1//2 ||F|| j [, 2 (c) 
follows from approximation of the solution in H 1 (O) by roughening. 

Similar results for the operator — div<7(x/e)V in a bounded domain with 
the Dirichlet or Neumann boundary condition were obtained in the papers 
[Grl, Gr2] by G. Griso by the "unfolding" method. 

0.3. Approximation of the resolvent in the {Li — > H 1 )-norm. The 

present paper relies on the results of [PSu]. In that paper, matrix DO's Ad, £ 
in a bounded domain O C M. d of class C 2 were studied. The operator Ad, £ 
is defined by the differential expression &(D)*g(x/e)6(D) with the Dirichlet 
condition on dO. The effective operator A® D is given by the expression 
b(D)*g°b(D) with the Dirichlet boundary condition. The behavior for small 
e of the solution u e of the equation Ad,s u s = F, where F G L2(C;C n ), is 
studied. Estimates for the i? 1 -norm of the difference of the solution u e and 
its first order approximation are obtained. By roughening of this result, an 
estimate for ||u £ — Uo||i, 2 (e>) is proved. Here uo is the solution of the equation 
A° D u = F. 
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In operator terms, the following estimates are obtained: 

\\A- D ) e - {A%r l - eK D {e)\\ L2{0) ^ m{0) < Ce 1 / 2 , (0.3) 

\\A D ] £ - (Alr l \\ L2{0) ^ L2{0) < Ce 1 ' 2 . (0.4) 

Here K £>(e) is the corresponding corrector. 

The method of [PSu] is based on using estimates (0.1), (0.2) for homoge- 
nization problem in M d obtained in [BSul,3] and on the tricks suggested in 
[Zh2], [ZhPas] that allow one to deduce estimate (0.3) from (0.1), (0.2). Main 
difficulties are related to estimating of the "discrepancy" w e , which satisfies 
the equation A £ w £ = in O and the boundary condition w £ = eKr>(e)F on 
dO. 

0.4. The main result. It must be mentioned that estimate (0.4) is quite 
a rough consequence of (0.3). So, the refinement of estimate (0.4) is a 
natural problem. In [ZhPas], for the case of the scalar elliptic operator 
— div<7(x/e)V (where g{x) is a matrix with real entries) an estimate for 

\\A~jJ~ £ - (A° d )~ 1 \\l 2 ^L2 °f order e 2d ~ 2 for d > 3 and of order e|loge| for 
d = 2 was obtained. The proof essentially relies on using the maximum 
principle which is specific for scalar elliptic equations. 

In the present paper, we prove a sharp order operator error estimate 

\\A~ D ) e - (^)- 1 ||l 2 ( OH l 2(0 ) < Ce. (0.5) 

Estimate (0.5) for matrix elliptic DO's refines even the known classical (non- 
operator) error estimates. 

Method of the proof relies on the results and technique of [PSu]. The 
problem reduces to estimating of the L2-norm of w e . Using of the operator 
approach and duality arguments is important. Employing approximation 
of the resolvent A~j D 1 £ in the norm of operators acting from L2(0;C n ) to 
Hq^O'jC 11 ), we find approximation of the same operator in the norm of op- 
erators acting from H~ l (0] C n ) to L2(0; C ra ). The last approximation com- 
bined with the boundary layer estimates allows one to obtain the required 
estimate for the I/2-norm of w e . 

0. 5. The plan of the paper. The paper contains three sections. In Section 

1, the class of operators is introduced, the effective operator is described, and 
the main result is formulated. Section 2 contains some auxiliary statements 
needed for further investigation. In Section 3, the main result is proved. 

0.6. Notation. Let fj and be complex separable Hilbert spaces. The 
symbols (•, and || • ||^ stand for the inner product and the norm in $); the 
symbol || • ||^^» denotes the norm of a linear continuous operator acting 
from f) to 

The symbols (•, •) and | • | stand for the inner product and the norm in C n ; 
1 = l n is the identity (nxn)-matrix. We use the notation x = (#i, . . . , x^) € 
R d , iDj = dj = d/dxj, j = l,...,d, D = -iV = (D 1 ,...,D d ). The L p - 
classes of C n -valued functions in a domain O cM. d are denoted by L p (0; C n ), 
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1 < p < oo. The Sobolev classes of C n -valued functions in a domain O C 
R d are denoted by H s (0;C n ). By H^(0;C n ) we denote the closure of 
Cg°(0;C n ) in i^OjC"). If n = 1, we write simply L p (0), H s (0), etc., 
but sometimes we use such abbreviated notation also for spaces of vector- 
valued or matrix- valued functions. 

0.7. Acknowledgement. The author is grateful to A. I. Nazarov for 
fruitful stimulating discussions. 

§1. Statement of the problem. Results 

1.1. The class of operators. Let Y C M. d be a lattice, and let O C l d be 

the elementary cell of the lattice Y. We denote |fi| = measfi. Below 
stands for the subspace of functions in H 1 ^) whose T-periodic extension to 
R d belongs to H^ oc (R d ). If p(x) is a T-periodic function in M d , we denote 

f £ {x) := <£(e _1 x), e > 0. 

Let O C R d be a bounded domain of class C 2 . In L2(0; C"), we define 
an operator Ad,z formally given by the differential expression 

Ae = 6(D)V(x)6(D) (1.1) 

with the Dirichlet condition on dO. Here g(x) is a measurable (m x m)- 
matrix- valued function (in general, with complex entries). We assume that 
g(x) is periodic with respect to the lattice Y, bounded and uniformly positive 
definite. Next, 6(D) = Ya=i b i D i 

is an (m x n)-matrix first order DO with 
constant coefficients. Here 6; are constant matrices (in general, with complex 
entries). The symbol b(£) = Yli=i h€h £ e corresponds to the operator 
6(D). It is assumed that m > n and that rank6(£) = n, V£ ^ 0. This 
condition is equivalent to the following inequalities 

a l n < b(0)*b(0) < ail„, 6 G < a < a x < oo, (1.2) 

with some positive constants «o and cci. 

The precise definition is the following: Ad,e is the selfadjoint operator in 
L2(0;C™) generated by the quadratic form 

a D , e [u,u] = / («f (x)6(D)u,6(D)u) dx, u G £#(0;C n ). 

Under the above assumptions this form is closed in L2(0;C n ) and positive 
definite. Moreover, we have 

c / |Du| 2 dx < a D;£ [u,u] < a / |Du| 2 dx, u G Hq(0; C"), (1.3) 

where Co = c^o Hi/ 1 llx^ 1 ' Cl = a illfl1Uoo- ^ i s eas Y to check (1.3) extending 
u by zero to M. d \ O, using the Fourier transformation and taking (1.2) into 
account. 

The simplest example of the operator (1.1) is the scalar elliptic operator 
Ae = — divg e (x)V = D*g £ (x)D. In this case we have n = 1, m = d, 
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6(D) = D. Obviously, condition (1.2) is valid with cto = cx\ = 1. Another 
example is the operator of elasticity theory which can be written in the form 
(1.1) with n = d, m = d(d + l)/2. These and other examples are discussed 
in [BSul] in detail. 

Our goal is to find approximation for small e for the operator Af ) 1 £ in 
the operator norm in L2{0;C n ). In terms of solutions, we are interested 
in the behavior of the generalized solution u e G i?Q(0;C n ) of the Dirichlet 
problem 

&(D)V(x)6(D)u e (x) = F(x), xeO; u £ j 9o = 0, (1.4) 
where F € L 2 (0; C n ). Then u e = A D ] £ F. 

1.2. The effective operator. In order to formulate the results, we need 
to introduce the effective operator A° D . 

Let an (n x m)-matrix-valued function A(x) be the (weak) T-periodic 
solution of the problem 

6(D)*«7(x) (6(D)A(x) + l m ) = 0, / A(x) dx = 0. (1.5) 

Jn 

In other words, for the columns Vj(x), j = 1, . . . , m, of the matrix A(x) the 
following is true: Vj € i? 1 ($7;C n ), we have 

/ ( 5 (x)(6(D)v,(x) + e J ),6(D)T 7 (x))dx = 0, Vij £ i? 1 ($7; C n ), 
Jn 

and J n Vj(x) dx. = 0. Here ei, . . . , e m is the standard orthonormal basis in 
C m . 

The so-called effective matrix g° of size mxmis defined as follows: 

g° = l^r 1 f 5 (x) (6(D)A(x) + l m ) dx. (1.6) 

It turns out that the matrix (1.6) is positive definite. The effective operator 
A° D for Ad,s is given by the differential expression 

-4° = 6(D)* 5 °6(D) 

with the Dirichlet condition on dO. The domain of this operator is 
Hl{0- C n ) n H 2 (0; C n ) (see Subsection 2.2 below). 
Consider the "homogenized" Dirichlet problem 

6(D)* 5 &(D)u (x) = F(x), x£0; u \ 9o = 0. (1.7) 

Then uo = (^2))" 1 F. As e —> 0, the solution u e of the problem (1.4) 
converges in Li{0\ C n ) to uo; for operators of the form (1.1) this was proved 
in [PSu]. We wish to estimate ||u e — Uo||l 2 (o). 

1.3. The main result. Denote 

{dO) £ = {xel d : dist {x, dO} < e}. 
Now we formulate the main result. 



6 T. A. SUSLINA 

Theorem 1.1. Assume that O C M. d is a bounded domain of class C 2 . 
Let g(x) and b(D) satisfy the assumptions of Subsection 1.1. Let u £ be the 
solution of the problem (1.4), and let uq be the solution of the problem (1.7) 
with F £ L 2 (0;C n ). Let e 1 £ (0,1] be such that the set {dO) £l can be 
covered by a finite number of open sets admitting diffeomorphisms of class 
C 2 rectifying the boundary dO. Let 2r\ = diamO, e 2 = + ri)~ 1 , and 
£0 = £2/2. Then for < e < eo we have 

ll u £ - Uo||l 2 (0;C") < Cis\\F\\ L , 2 ( . C n), (1.8) 

or, in operator terms, 

\\^D,e ~ {Ad)~ 1 \\l 2 {0;C™)^L 2 {0;C") < C\£. 

The constant C\ depends only on m, d, ao, a±, H^Hl^, ||<7 _1 ||.Loo! the pa- 
rameters of the lattice T, and the domain O. 

§2. Auxiliary statements 

2.1. The energy inequality. Consider the problem (1.4) with the right- 
hand side of class H- l {0;C n ). Recall that H- l {0;C n ) is defined as the 
space dual to H^(0;C n ) with respect to the L 2 (0; C n )-coupling. If f £ 
fr _1 (0;C n ) and 77 £ H%(0;C n ), then the symbol ({,v)l 2 (0) = J (i,r))dx 
stands for the value of the functional f on the element rj. Herewith, 



f , rj) dx. 

o 



< llfl 



H- 1 (0;C n )\\ r l\\H 1 (0;C n )- 



The following (standard) statement was checked in [PSu, Lemma 4.1]. 
Lemma 2.1. Let f £ H~ 1 (0;C n ). Suppose that z e £ H^(0;C n ) is the 
generalized solution of the Dirichlet problem 

6(D)V(x)6(D)z £ (x) = f(x), x£0; z £ \ do = 0. 

Ln other words, z £ satisfies the identity 

[ (<f(x)6(D)z £ ,6(D)r;)dx = / (f,r,)dx, V r, £ H^O; C n ). 
Jo Jo 

Then the following "energy inequality" is true: 

\\ z e\\m(0;C n ) < C||f ||fl--i(o ; C")) 
where C = (1 + (diam O) 2 )a 1 \\g~ l \\ Loo . 

It follows from Lemma 2.1 that the operator A~ r ] 1 £ acting in L 2 (0;C n ) 
can be extended to a linear continuous operator acting from i7 _1 (C; C n ) to 
i7g(C;C n ). Applying Lemma 2.1 with g 6 replaced by g°, we see that the 
same statement is true for the operator (A^)^ 1 . 

Note that 

(A D ) £ h,h) L2{ o) = (fuA^^o), fi,f 2 G H~\0-C n ). (2.1) 
A similar identity is valid for the operator (A^)^ 1 . 
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All the statements of Subsection 2.1 are valid in arbitrary bounded domain 
O (without assumption that dO G C 2 ). 

2.2. Properties of the solution of the homogenized problem. Due 

to the assumption dO G C 2 , the solution uo of the problem (1.7) satisfies 
u G Hl{0-C n ) DH 2 (0;C n ), and 

ll u \\h 2 (0; C") < c|| F IIl 2 (C;C")- ( 2 - 2 ) 
In operator terms, it means that 

ll(>^£)) _1 |lL2(0;C«)^if 2 (0;C™) <c. (2.3) 

The constant c depends only on cto, a±, H^Hl^ , Hfl 1-1 Hi^ , and the domain O. 
To justify these properties, it suffices to note that the operator b(D)*g°b(D) 
is a strongly elliptic matrix DO and to apply the "additional smoothness" 
theorems for solutions of strongly elliptic systems (see, e. g., [McL, Chapter 
4])- 

2.3. Trace lemma. We need the following simple statement; see, e. g., 
[PSu, Lemma 5.1]. 

Lemma 2.2. Denote B £ = {x G O : dist {x, dO} < e}. Then for any z G 
H\0) we have 

I \z\ 2 dx < (3e\\z\\ H i( 0) ||z||l 2 (o), < e < e ± . 
Jb £ 

Here £\ is the same as in Theorem 1.1. The constant (3 depends only on the 
domain O. 

Note that the statement of Lemma 2.2 is valid for any bounded domain 
Oof class C 1 . 

2.4. Smoothing in Steklov's sense. Let S £ be the operator in 
L 2 (R d ;C m ) given by 

(S e u)(x) = \n\- 1 [ u(x-ez)dz. (2.4) 
Jn 

It is said that the operator S £ is smoothing in Steklov's sense. 

We need the following property of the operator (2.4) (see [ZhPas, Lemma 
1.1] or [PSu, Proposition 3.2]). 

Lemma 2.3. Let /(x) be a T-periodic function in M. d such that f G ^(O). 
Let [f £ ] denote the operator of multiplication by the function / £ (x). Then 
the operator [f £ ]S £ is continuous in L2(M d ;C m ), and 

\\[f £ ]Se\\L2(R d ;C™)^L 2 (m. d ;C™) < 1^1 1/2 ll/llx, 2 (n) • 

2.5. Properties of the matrix A(x). Let T be the lattice dual to T. By 
O we denote the central Brillouin zone of T, i. e., Q = {k G M. d : |k| < 
|k — b|, / b G T}. Let ro be the radius of the ball inscribed in clos Q. 
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Recall that the matrix-valued function A(x) is the T-periodic solution of 
the problem (1.5). In [BSu2, Subsection 7.3] it was proved that 

II a II ^ 1/2 /o \-llo|l/2 — 1 /2 1| n 1/2 — I,, 1/2 (n _s 

||A|| La (n) < m ' (2r ) 7 a IbH^Jb llz^- ( 2 - 5 ) 

Let [A e ] be the operator of multiplication by the matrix-valued function 
A £ (x); this operator acts from L 2 (R d ;C m ) to L 2 (R d ;C n ). By Lemma 2.3 
and estimate (2.5), the norm of the operator [A e ]S £ satisfies the following 
estimate: 

\\[h £ }Ss\\L 2 (M. d ;C m )^L 2 (R d ;C n ) ^ 1^1 ^ II M\l 2 (Q) 
s l/2/ n-1 -1/2m ||1/2|| —1 1| 1/2 „, ^ - b ) 

<m> (2r ) a Hsll/Jb IIl 00 = :M - 

§3. Proof of Theorem 1.1 

The proof of Theorem 1.1 relies on the results of [PSu], where approxima- 
tion of A~j^ £ in the norm of operators acting from L 2 (0;C n ) to H l (0;C n ) 
was obtained. 

3.1. Error estimates in H 1 . We fix a linear continuous extension operator 
P a ■ H 2 (0; C n ) -> H 2 (R d ; C n ) (3.1) 
and put uo = Po u o- Then 

ll u oll_H' 2 (R d ;C n ) — Co\\ u o\\H 2 (0;C n )' (3-2) 

where Co is the norm of the operator (3.1). Let S e be the smoothing 
operator (2.4). By Rq we denote the operator of restriction of functions in 
M. d onto the domain O. We put 

K D {e) = RoW]S e b{T>)Po{A Q D )-\ (3.3) 

The operator b(D)Po{A D )~ 1 is a continuous mapping of L 2 (0;C n ) into 
H 1 (R d ;C m ). Using Lemma 2.3 and relation A € H 1 ^), it is easy to check 
that the operator [A £ ]S £ is continuous from il^R^C™) to ff 1 ^; C n ). 
Hence, the operator (3.3) is continuous from L 2 (0;C n ) to H l {0]C a ). 
The following statement was proved in [PSu, (7.10)]. 

Proposition 3.1. Let O cM. d be a bounded domain of class C 2 . Let u £ be 
the solution of the problem (1.4), and let uq be the solution of the problem 
(1.7) with F G L 2 (0;C n ). Xei uo = -Pe>uo, where Po is the extension oper- 
ator (3.1). Let w £ G H l {0]C n ) be the generalized solution of the problem 

Aw £ = 0inO, w £ | 90 = eA £ S £ &(D)uo|do- (3.4) 

T/ien for < e < 1 we have 

||u £ - u - eA e S , e 6(D)u + w £ || H i (0 . c „) < CellF^^.cn). (3.5) 

The constant C depends only on m, d, ao, a±, H^Hloo, ||5 _1 ||l 00) the param- 
eters of the lattice T and the domain O. 

The following theorem was proved in [PSu, Theorem 7.1]. 



HOMOGENIZATION OF THE ELLIPTIC DIRICHLET PROBLEM 9 

Theorem 3.2. Suppose that the assumptions of Theorem 1.1 are satis- 
fied. Let uo = -Pe>uo, where Pq is the extension operator (3.1). Then for 
< e < £ 2 we have 

||u £ - u - eA £ S £ b(U)u \\ H i {o . Cn) < Ce^llFHi^c-cn), (3.6) 

or, in operator terms, 

W^D,e ~ (Ad)' 1 ~ £K D (£)\\ L2 ( . C n^ H i( . Cn) < Ce 1/2 . 

The constant C depends only on m, d, ceo, ct±, \\gl\Loo, ||<7 _1 ||.Loo! ^ e param- 
eters of the lattice T, and the domain O. 

Recall that {dO) £ denotes the e-neighborhood of dO. For sufficiently 

small e, we fix two cut-off functions £ (x) and £ (x) in W 1 such that 

9 £ eC%°{R d ), supple (dO) £ , O<0 £ (x)<l, 
# £ ( x )|de> = 1) £ | V6' e (x) | < k = const; 



9 £ eC^°(M. d ), supp^C (30) 2£ , O<0 £ (x)<l, 



< k = const. 



(3.7) 
(3.8) 



£ (x) = 1 for xe (dO) £ , e V# e (x) 
We denote 

<f> £ = £6 £ A £ S £ b(T>)u . (3.9) 
From (1.2), (2.2), (2.6), (3.2), and (3.7) it follows that 

II0 £ IU 2 (O;C™) < *i ll u o||Hi(Rd;C") < £Ma i Ce»cll F IU 2 (0;C")> ( 3 - 10 ) 

cf. [PSu, (7.14)]. The norm of the function (3.9) in H l (0; C n ) was estimated 
in [PSu, Lemma 7.4]. A similar estimate is true if 6 £ is replaced by 9 £ . We 
formulate the corresponding result. 

Lemma 3.3. Suppose that the assumptions of Theorem 1.1 are satisfied. 
Let 9 £ and 9 £ be functions satisfying (3.7), (3.8), and let <f> £ be defined by 
(3.9). Then we have 

\\4>e\\H^{0;C n ) < C 2 £ 1/2 ||F|| i2 ( 0;C n) , < £ < £ 2 , (3.11) 

\\£e £ A £ S £ b(B)u \\ H1{o . Cn) < C 2 £ 1 / 2 \\F\\ L2{0 . Cn) , < 2e < e 2 . (3.12) 

The constants C 2 and C 2 depend only on m, d, ao, ai, HsHi^, ||5 _1 ||loo5 
the parameters of the lattice T, and the domain O. 

3.2. Proof of Theorem 1.1. Step 1. Roughening (3.5), we obtain 

||u £ - u - eA e S e 6(D)uo + w £ || L2(0;C n) < Ce\\F\\ L ^ . C n), < e < 1. 

(3.13) 

Combining (1.2), (2.2), (2.6), and (3.2), we see that 

||A £ S £ 6(D)S ||l 2 ( O ) < Ma^CocWnL^oy (3-14) 
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From (3.13) and (3.14) it follows that 

K " u oIIl 2 (o) < e(C + Ma\ /2 C cW\\L 2 (0) + l|w £ || MO) , < e < 1. 

(3.15) 

Therefore, the proof of estimate (1.8) is reduced to estimating of w £ in 
L 2 (0;C n ). 

For this purpose, we need the following lemma. 
Lemma 3.4. Suppose that the assumptions of Theorem 1.1 are satisfied. 
Let 9 £ be a function satisfying (3.8). Consider the operator 

K D {e) = R [(l - 9 £ )A £ }S £ b(D)P o (A D )-\ (3.16) 

which is a continuous mapping of L 2 (0;C n ) into -£^(0; C ra ). Let (Kr>{e))* : 
if~ 1 (0;C n ) L 2 (0;C n ) be the operator adjoint to the operator (3.16), i. e., 

V Jl 2 (o) V )l 2 {o) (3. 17 ) 

Vf €fl" _1 (0;C n ), VveL 2 (C;C n ). 

Then the operator A : D 1 £ , viewed as a continuous mapping of .ff _1 ((9; C n ) 
into L 2 (0;C n ), admits the following approximation 

\\A D \ - {Air 1 - e(K D {e)r\\ H - HO) ^ L2io) < (C + C 2 )e 1 '\ < 2e < e 2 . 

(3.18) 

Proof. From (3.6) and (3.12) it follows that 

K - u - e(l - e £ )A £ S £ b(-D)u \\ HHo . cn) <(C + C 2 )e 1 / 2 ||F|| L2(0 . C „ ) , 

< 2e < e 2 . 

_ (3-19) 

The function under the norm-sign on the left belongs to H^(0;C n ). In 
operator terms, (3.19) means that 

WA^-iAD^-ekDie^o^gi^ < (C+C^e 1 / 2 , 0<2e<e 2 . (3.20) 

This implies (3.18) by the duality arguments. Indeed, combining (2.1), the 
similar identity for (A^^ 1 and (3.17), we see that for any f G i? _1 (C;C n ) 
and v G L 2 (0; C n ) one has 



((A D ] £ - (Al)' 1 - e(K D (e)r)i,v) 
= (f,(A D ] £ -(A° D r 1 -ek D (e))v) 



L 2 (0) 



L 2 (0) 

Together with (3.20) this yields 

((^; £ -(^)- i - e (^( £ ))*)f,v) L2(o) 

< (C + C 2 )e 1/2 ||f \\ H - HO) ||v|| i2(0) , Vf €fr 1 (0;C n ), V v € L 2 (0; C n ), 
which implies (3.18). • 
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From (3.7) and (3.9) it follows that 4> £ \eo = eA E <S e &(D)uo|a0. Then, by 
(3.4), the function w £ — <f> £ is the solution of the problem 

A e (w e - <f> e ) = F £ in O, (w £ -0 £ )| 9o = 0, (3.21) 

where F £ = —A £ 4> £ . It is easily seen that F £ G //~ 1 (C;C n ), and 

HFellif-HOjC™) < a l rfl/2 !lfi'l|Looll ( / , ell J H' 1 (C';C")) (3.22) 

see [PSu, (4.15)]. From (3.11) and (3.22) it follows that 

I|F 6 ||h-i(0;C») < C3e 1/2 ||F|| L2(0;C n), < e < e 2 , (3.23) 

where C3 = a\d l / 2 \\g \\l oc C 2 . Note also that F £ is supported in (dO) £ . 

Now we apply approximation (3.18) to the problem (3.21). Since w £ — 
(f> £ = A^ £ F £ , then 

||w £ - 4>e ~ (Ar'Fe - e(K D (e))*F £ \\ L2{0) < (C + C 2 )e 1/2 ||F e || H -i (0)) 

< 2e < e 2 . 

(3.24) 

By (3.16) and (3.17), for any v G L 2 (0;C n ) we have 
((K D (e))*F £ ,v) ^ = (F £ ,(l-9 £ )^S £ b(B)P o (A D )- 1 v) L ^ oy (3.25) 

Since 1 — £ (x) = for dist {x, dO} < e, and F £ is supported in (dO) £ , then 
the right-hand side of (3.25) is equal to zero. Consequently, (Kd(e))*F £ = 0. 
Then (3.24) and (3.23) imply that 

||w £ -<f> £ - (Alr^shzio) < (C7 + C , 2 )C 3 e||F|| L2(0) , < 2e < e 2 . (3.26) 

The norm of cb £ in L 2 (0; C n ) admits estimate (3.10). It remains to estimate 
the L2-norm of the function (A ( £ ) )~ 1 F £ . 

3.3. Proof of Theorem 1.1. Step 2. 

Lemma 3.5. Suppose that the assumptions of Theorem 1.1 are satisfied. 
Let 4> £ be defined by (3.9), and let F £ = -A £ (j> £ = -6(D)* g £ b(T))(f> £ . Then 
the function rj £ := (^4j)) _1 F £ satisfies the following estimate: 

\\Ve\\L 2 (0;C n ) < C^HFU^o-cn), < £ < £ 2 . (3.27) 

The constant C4 depends only on m, d, Qo, «i, HpIIlooj II<7 _1 ||looi the pa- 
rameters of the lattice T, and the domain O. 

Proof. The function rj £ € ^(OjC™) is the generalized solution of the 
Dirichlet problem A°r/ £ = F £ , r/ £ \Qo = 0. It means that 

/ ( 5 °6(Dh £ ,6(D)h)dx= / (F £ ,h)dx = - / (/6(D)0 £ ,6(D)h)dx, 
Jo Jo Jo 

Vhe Hl{0-<C n ). 

(3.28) 
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If h € H 2 (0;C n ) fl H^(0;C n ), then it is possible to integrate by parts in 
the left-hand side of (3.28). Hence, 

I (rj £ ,A°h)dx=- f (/6(D)0 £ ,6(D)h)dx, 
Jo Jo (3.29) 

Vh G H 2 (0;C n ) f\Hlip-C n ). 

Now we write down the norm of the function 7] £ in -^(C; C n ) as the norm 
of continuous antilinear functional: 

|/ (77 £ ,G)dx| 

Nellz^o) = SU P — mFhi • 

0^GeL 2 (O;C") IMIl 2 (0) 

We put h = (A^^G, G G L 2 (C; C n ). Then G = A°h, and h runs through 
H 2 (0;C n ) n H^(0;C n ) if G runs through L 2 {0;C n ) (see Subsection 2.2). 
Hence, 

\Uri e ,A h)dx.\ , s 

\\Ve\\L 2 (0)= SUp • (3 - 30) 

0^hefl-2(O)nHi(O) II- 4 n llL 2 (o) 

By (2.3), we have ||-A h|| L2(o) > (cT 1 ||h|| H 2 (cl) . Combining this with (3.29) 
and (3.30), we obtain 

I, I, ^ |/ o (^ft(D)0 e ,&(D)h)rfx| 

\\Vs\\l 2 (o) < c sup ^ — -. (3.31) 

o^hetf 2 (O)n^(O) ll n ll// 2 (e>) 

Next, since 6(D) = Y,i=i b i D i and > b Y ( L2 )> N < «i /2 > then 

||/6(D)0 e || L2(o) < ||5||Loc«i /2 ^ 1/2 ||0 £ lki(O)- (3-32) 

Taking into account that the function <j> £ is supported in the e-neighborhood 
of dO, from (3.31) and (3.32) we see that 

- ~. „ 1/2 ,1/2,, , „ ll & ( D ) h lli«(B e ) 

\\Ve\\L 2 (0) < C||£f|koo«l O 7 ||0 e ||Hi(O) SUp — . 

0^he-ff 2 (O)n^i(O) ll n !lff 2 (e>) 

(3.33) 

Applying Lemma 2.2 and taking into account that |6(D)h| < 
a 1 / X)j=i |Ah|, for < e < ei we have: 



/ \b(D)h\ 2 dx < ai d^T f \Dih\ 2 dx 

J B e JB £ 



< o;id/3e^ ||Ah||Hi(o)||Ah||L 2 (o) < aid/3e||h|| H 2 (cl) ||h|| H i (C)) . 
1=1 



Hence, 



l|6(D)h|| L2(Be ) < . , M o/ 2 ,/ 2 



sup " < (aiWe 1 ^, < e < £1. (3.34) 

0^heif 2 (O)nifi(O) ll h lltf 2 (e>) 
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Finally, from (3.11), (3.33), and (3.34) it follows that (3.27) is valid with 
C4 = c\\g\\ Lao P 1/2 a 1 dC 2 . • 

Now it is easy to complete the proof of Theorem 1.1. By (3.10), (3.26), 
and (3.27), we have 

II w e ||l 2 (0) < Cke||F|U 2 (o), < 2e < e 2 , 

where C 5 = (C + C 2 )C 3 + Ma] /2 C c + C 4 . Combining this with (3.15), we 
arrive at (1.8) with d = C + Ma\ /2 C c+ C 5 . • 
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